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Abstract: Lattice data, taken since 1998 near the critical line of a 4D Ising model, 
have been supporting the large logarithmic rescahng of the scalar condensate predicted in 
the alternative description of symmetry breaking proposed by Consoli and Stevenson. This 
conclusion has been challenged in a recent paper by Balog et al. In this paper we respond to 
the criticism of these authors, recapitulate the theoretical and numerical evidences in favour 
of the alternative interpretation of 'triviality' and reiterate our conclusion: 'triviality', by 
itself, cannot be used to place upper bounds on the Higgs boson mass. 
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1. Introduction 

Recent lattice data [1], collected near the critical line of a 4D Ising model, support the 
large logarithmic rescaling of the scalar condensate predicted in an alternative description 
of symmetry breaking in theories, see Refs. [2-4]. This result, while confirming previous 
numerical indications obtained in Refs. [5-7], would have a substantial phenomenological 
implication: one cannot use 'triviality' to place upper bounds on the Higgs boson mass. 

This point of view has been challenged in a recent paper [8] by Balog et al.. These au- 
thors, referring just to Ref. [1], while otherwise ignoring the previous numerical indications 
of Refs. [5-7], draw the opposite conclusion: the standard interpretation, as they say the 
'Conventional Wisdom' (CW), is completely consistent with all lattice data. The aim of 
this paper is to respond to their criticism, recapitulate in a unified framework the results 
of Ref. [1] and Refs. [5-7], and reiterate our conclusion: 'triviality', by itself, cannot be 
used to place upper bounds on the Higgs boson mass. 

2. The rescaling of the scalar condensate 

Before entering the details of the controversy, we shall first remind once again in this section 
why in a spontaneously broken phase there are two basically different definitions of the field 
rescaling. In fact, the widespread skepticism concerning the interpretation of 'triviality' 
proposed in Refs. [2-4] originates from a non sufficient appreciation of this crucial point. 

To this end, let us introduce the bare 'lattice' field ^b{x) = ^\a,t\.{x) (i.e. as defined 
at a locality scale fixed by the ultraviolet cutoff A ~ 7r/a, o being the lattice spacing) and 
its expectation value (the 'scalar condensate') 
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Connecting to the stability analysis, the values ±vb represent the absolute minima of 
the effective potential Veff{^pB) of the theory. We shall also introduce the bare shifted 
fluctuation field 

hB{x) = ^B{x)-{^B), (2.2) 

whose expectation value vanishes by definition. 

Now, a first natural definition of the field rescaling, say Z = Zprop, is obtained from 
the residue of the shifted-field propagator 

G.oUs) - (2.3) 

near the physical mass-shell s = m\ . This can be used to define a renormalized fluctuation 
field hR{x) 

hB{x) = \/Zp^ophR{x) (2.4) 

whose propagator has the canonical form. 

This first definition of the field rescaling is constrained by the Kallen-Lehmann repre- 
sentation to lie in the range < Zp^op < 1, the free-field limit corresponding to the case 
^pvop = 1- This can be rigorously established for the local, Lorcntz-covariant theory. If 
this is viewed as the continuum limit of the cutoff theory, one expects [9] 

Z^rlp = 1 + / dspAis) (2.5) 

J So 

where Sq denotes the continuum threshold (sq = C^f^H)'^ in perturbation theory) and pa{s) 
the spectral function. For this reason, in the continuum limit A ^ oo, where according to 
'triviality' the spectral function should tend to 6{s — mjj), Zp^op should tend to unity. 

Another definition of the field rescaling, say Z = Z^, is peculiar of a broken-symmetry 
phase. It indicates the rescaling that is needed to relate the physical vacuum field vr to 
the bare vb, i.e. 

= -7^ ■ (2-6) 

By physical, we mean that the second derivative of the effective potential V^f^^ipR) evaluated 
at ifR = ivR, is precisely given by the physical Higgs boson mass squared, i.e. 

Kff(^K) = ml (2.7) 

This is very simple to understand. V^q{(Pr) represents the renormalized 2-point function 
evaluated at an external 4-momentum = 0. For ipR = vr, this should match the inverse 
of the renormalized connected propagator at zero momentum. In a 'trivial' theory, in the 
continuum limit, this has the simple free-field form Gr{0) = l/mj^. 
Therefore, this other definition is equivalent to the relation 



K4(-b) 



Z^ = T7FT:^=mlx2{0) (2.8) 



where X2(0) = 1/V^g(^s) is the bare zero-momentum susceptibility. 
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Let us now explain why Zprop and Z^p are completely different physical quantities. To 
this end, we shall consider the class of 'triviality compatible', gaussian-like approximations 
to the effective potential, say T4fr('/'_B) = V^trivCv^s)) where the shifted fluctuating field is 
governed by an effective quadratic hamiltonian. In this class, either Zprop = 1 identically 
or it tends to unity in the continuum limit A — > cxd. In spite of this, the rescaling of the 
vacuum field 

^^ = "i|/X2(0)~lnA (2.9) 

diverges logarithmically. 

This result is easily recovered noticing that the class Vtriv includes the one-loop poten- 
tial, the gaussian approximation and the infinite set of post-gaussian calculations where 
the effective potential reduces to the sum of a classical background energy and of the zero- 
point energy of the free massive shifted field with a c/^^— dependent mass, ^2 = i7((/?^). For 

instance, at one loop one finds V?{ip\) = m^-|-— and \b denoting respectively the 
bare mass squared and the bare fourth-order self-coupling. In gaussian and post-gaussian 
calculations, the functional dependence = J7((^^) is determined upon a minimization 
procedure of the energy functional in a suitable class of quantum trial states. In all cases, 
the basic mass parameter of the broken-symmetry phase is obtained through the relation 
ran = ^(vb). 

Therefore, the various approximations belonging to Vtriv share the same simple general 
structure, up to non-leading terms that vanish faster when A — > oo (see Ref. [4] for the 
general argument, Ref. [10] for the gaussian approximation and Ref. [11] for the case of 
a post-gaussian approximation). For the particularly simple case of the classically scale- 
invariant theory (the 'Coleman- Weinberg regime') this is given by 

Aeff^l, ^]^(^|) C7^(^|) 1 

VtriviVB) = —rr- + To. n) (2-10) 



4! 647r2 ' A2 2' 



with 



^H^b) = ^ (2.11) 

all differences among the various approximations being isolated in the relation between the 
effective coupling Agg = Aefr(O) and the bare coupling. 

For instance, at one loop Xes = Xb while in the gaussian approximation (see Ref. [10]) 
Aefif corresponds to resum all one- loop bubbles with mass = ^{(p%) so that 



Now, since all approximations display the same general structure in Eqs.( |2.10| ) and (2.11), 



upon minimization and after setting (fB = ^vb, one finds the same leading logarithmic 
trend of the effective coupling 

Xcsimn) ~ -— r (2-13) 
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Formally, this is precisely the same relation obtained in perturbation theory from the 
leading-order (LO) renormalized coupling 

after sending ^ oo at the Landau pole 

^ . ^ 16vr2 

Koimu) ~ (2.15) 
31n 

Therefore, one might be tempted to conclude that, at least at the leading logarithmic level, 
the alternative picture of Refs. [2-4], is equivalent to the CW. 

However, there is one notable difference: the quadratic shape of the 'trivial' effective 
potential Eq.( 2.10 ) at its absolute minima, namely the inverse susceptibility V^riv(^s), is 
not given by = 17^ (f^). Rather, one finds 

2 

<iv(^B) ~ Aefrml, ~ (2.16) 

thus leading to Eq.(^). 

Eq. ( p.g| ) has substantial phenomenological implications. In fact, using Eq. ( |2.6D with 
a logarithmically divergent Z^p as in Eq. ( |2.9| ), one finds 

(2.17) 

instead of the conventional relation obtained for a trivial unit rescaling 

[<]c^-^l (2-18) 
Now, whatever the vr — vb relation, the Higgs boson mass squared goes like 

ml, ~ (2.19) 



This becomes 



2 ["rI CW 10 on\ 
^^^J- (2.20) 



by using Eq.( |2.18D but becomes 

ml ~ 4 (2.21) 

if one defines vr through Eq.( p.l7 ). In the latter case, uih and vr scale uniformly in the 
continuum limit. 

Therefore, assuming to relate the vr of Eq.( |2.17] ) (and not the conventional [i^rJc^ ~ 
Vb oi Eq.( 2.18] )) to some physical scale (say 246 GeV), a measurement of uih would not 
provide any information on the magnitude of A since, according to Eq.( 2.21 ), the ratio 
C = tuh /vr is now a cutoff-independent quantity. Moreover, in this approach, the quantity 
C does not represent the measure of any observable interaction (see the Conclusions of 
Ref. [12]). 
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We emphasize that the difference between and Zprop has a precise physical mean- 
ing being a distinctive feature of the Bose condensation phenomenon [4]. In the class of 
'triviality-compatible' approximations to the effective potential, one finds C = mH/vji = 
27rv^2(', with < < 2 [A], ( being a cutoff-independent number determined by the 
quadratic shape of the effective potential Vesi^R) at = 0. For instance, C = 1 corre- 
sponds to the classically scale-invariant case or 'Coleman- Weinberg regime'. 

As for the standard interpretation of 'triviality', the gaussian effective potential ap- 
proach can also be extended to any number N of scalar field components [13] . In particular, 
when studying the continuum limit in the large-N limit of the theory, one has to take into 
account the non-uniformity of the two limits, cutoff A ^ oo and N ^ oo [11, 14]. This is 
crucial to understand the difference with respect to the standard large-N analysis. 



3. Lattice tests of the alternative interpretation of 'triviality'. Part 1. 

To test the alternative picture of 'triviality' of Refs. [2-4] against the CW one can run 
numerical simulations of the theory and check the scaling properties of the squared Higgs 
lattice mass against those of the inverse zero-momentum lattice susceptibility. 

The numerical simulations of Ref. [1] (and of Refs. [5-7]) were performed in the Ising 
limit that traditionally has been chosen as a convenient laboratory for the numerical anal- 
ysis of the theory. 

In this limit, a one-component $| theory becomes governed by the lattice action 

•Sising = -'^ X] X] ^'^^^ + e/,)(/)(a;) + 0(x - e^)(/){x)] (3.1) 

X fj, 

where (j){x) takes only the values ±1. 

Using the Swendsen-Wang [15], and Wolff [16] cluster algorithms we computed the 
bare magnetization: 

VB = M) , ct>^Y.^{x)/L^ (3.2) 

X 

(where (j) is the average field for each lattice configuration) and the bare zero-momentum 
susceptibility: 



Xlatt = L 



mr ■ (3.3) 



We report in Table 1 our determinations of vb and Xiatt- Other values, obtained over the 
years by different authors, are reported in Table 1 of Ref. [8]. 

Checks of the logarithmic trend predicted in Eq.( |2.9D can be performed using different 
methods. In a first indirect approach, discussed in this Section, one can use the fact that, 
both in perturbation theory and according to Refs. [2-4], the bare field expectation value 
is predicted to diverge logarithmically in units of the physical Higgs boson mass mn, 

i.e. 

~|ln(K-K,)|. (3.4) 
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Therefore, following CW, where the zero-momentum susceptibility is predicted to scale 
uniformly with the inverse squared Higgs mass 



hlfX2(0)]cw~l (3-5) 

one expects 

[4x2(0)]cw~IM'^-'^c)| (3.6) 

On the other hand, in the approach of Refs. [2-4] one predicts Z^p = X2(0)w.|^ ^ ln(A) so 
that, in this case, one would rather expect (CS=Consoli-Stevenson) 

[4x2(0)]cs~|ln('«-MP (3.7) 



The two leading-order predictions in Eq. ( p.6| ) and in Eq. ( p.T] ) were directly compared 
with the lattice data for the product f^X2(0) reported in Table 3 of Ref. [1]. These data 
were fitted to a 3-parameter form 

a|ln(K-Kc)r (3.8) 
where a is a normalization constant and one can set the exponent 7 = 1, according to 



Eq. (|3.6|) , or 7 = 2 according to Eq. (|3.7D . Using this type of functional form the results of 
the fit to the lattice data single out unambiguously the alternative picture of Refs. [2-4], 
i.e. 7 = 2, over the value 7 = 1 (see Fig. 2 of Ref. [1]). 

However, it has been pointed out by the authors of Ref. [8] that the lattice data can 
also be reproduced for 7 = 1, once the scale within the log is left as a free parameter, 
thus effectively simulating the presence of next-to-leading corrections. Indeed, in this case 
one gets a good fit in both cases with precise determinations of the critical point, namely 
Kc = 0.074833(17) for 7 = 1 and Kc = 0.074819(15) for 7 = 2, in good agreement with the 
value Kc = 0.074834(15) obtained by Gaunt et al. [17] from the symmetric phase. 

In this sense, we can agree with Balog et al.: a definitive test to decide between the 
two asymptotic behaviours 7 = 1 and 7 = 2 has to be postponed to simulations performed 
closer to the critical point where the non-leading terms associated with the scale of the 
logs should become unessential. 

However, although one can certainly get a good fit with 7 = 1, the agreement between 
lattice data and the prediction based on 2-loop renormalized perturbation theory is not 
good. This can be checked considering the expression reported by Balog et al. (/ = 

1 \Vl{k - Kc)|) 

25 

fe2(0)]2_i„„p = ai(Z - - In/) + a2 (3.9) 
together with the theoretical relations [8] 

(3.10, 

and 

— =ln(C^) + 2 ln(C;) -1.6317 (3.11) 
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Figure 1: We show the lattice data for w^xiatt of Table 1 together with the fit Eq. (3.8) for 7 = 2 



(solid line) and the two-loop fit Eq. (^^) where the fit parameters are allowed to vary inside their 
theoretical uncertainties, ai = 1.20(3) and 02 = —1.6(5) (dashed line). 



Using the input values reported by Liischer and Weisz (LW) in Table 1 of Ref . [18] and 
the relations C[ = e^'^Ci, C'^ = C2 and = C3 (see Eqs.(4.37) and (4.38) of Ref. [18]), 
the predictions for the Ising model (A = 1) are: = 6.49(7), 21n(C() = 1/3 + 3.0(4) and 
ln(C3) = —3.0(1) or ai = 1.20(3), 02 = —1.6(5). However, the data for f^X2(0) require 
oi = 1.267(14) and 02 = —2.89(8), see Ref. [8]. Therefore, the quality of the 2-loop fit is 
poor (see Fig.l). 

Nevertheless, one can adopt a pragmatic point of view, ignoring possible problems 
related to the matching conditions with the symmetric phase, and try to extract from the 
data for f^X2(0) a new set of constants C'^. In particular the precise determination from 
ai = 1.267(14) implies 

C2 = 6.67(4) (3.12) 

that will be used in the following. 

Thus we can summarize the results of this section as follows: 
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1) leaving out the scale within the logs as a free parameter (and thus allowing 
effectively for the presence of next-to- leading corrections) the lattice data for w^X2(0) are 
unable to distinguish between the powers 7 = 1 and 7 = 2. In this sense, a definitive 
test has to be postponed to simulations performed closer to the critical point where the 
non-leading terms associated with the scale of the logs should become unessential. 

2) a 2-loop fit as in Eq.(^^), with ai and a2 deduced consistently from the 
LW tables, does not provide a good description of the data (see Fig.l). However, within 
the effective 2-loop formula Eq.(^.S|), when ai and 02 are left as free parameters, one can 
obtain precise determinations of the critical point Kc = 0.074833(17) and of the integration 
constant C2 = 6.67(4). These will be used in the next section to check whether the lattice 
observables are consistent with the critical behaviour predicted by Renormalization Group 
(RG) analysis. 

4. Lattice tests of the alternative interpretation of 'triviality'. Part 2. 

Additional numerical evidences concerning the relative scaling of niH and X2(0) can be 
obtained by comparing again with the predictions of perturbation theory. To this end, 
we performed in Ref. [1] a test of the logarithmic trend predicted in Eq.( |2.9D assuming as 
input entries the full 3-loop values minput = mR reported in the first column of Table 3 of 
Ref. [18] at the various values of k. These input mass values, to leading order, follow the 
scaling law 

ruH = A^K - Kc-\ ln(fi; - Kc)\~^^^ (4.1) 
so that one can check whether the quantity 

= 2Kmfj,p^tXiatt (4.2) 

tends to unity or grows logarithmically when approaching the continuum limit. 

Now, using in Eq.( |4.2| ) the central values of miaput = reported in the LW Table 
and the values of xiatt reported in Table 1, the conclusion is unambiguous: the in 
Eq.([4.2|) becomes larger and larger approaching the continuum limit along the RG curve 
"T-input = fnmp\it{i^) and the observed increase is completely consistent with the logarithmic 
trend predicted in Eq.( |2.9|) (see Fig.2). 

The discrepancy with the perturbative predictions can also be checked noticing that 
for minput = 0.08 Ref. [18] predicts k = 0.07481(8), i.e. smaller than 0.0749. Therefore, 
by inspection of Table 1, the relevant lattice susceptibility will be definitely larger than 
its value for k = 0.0749, xiatt ~ 1100, so that using Eq. (^), one gets the lower bound 
> 1.05 which cannot be reconciled with the perturbative predictions. 

Balog et al. object to our conclusions that "..the crucial question is whether the 
estimates minput of vtlr are reliable". According to these authors, "..the measured values 
of rriR are considerably lower than the corresponding estimates minput" (i-e. the m/j's that 
we took from the LW table). As a matter of fact, by replacing the LW m/j's with the 
results of their simulations, one gets a remarkably constant value 2Km|jXiatt ~ 0.88. 

Thus their objection does not concern our strategy but rather the validity of the LW 
entries themselves as reliable estimates of the Higgs mass parameter m//. Actually, as we 
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^input 


K 


lattice 


algorithm 


Ksweeps 


Xlatt 


/ 1 J. 1 \ 


n A 

OA 


0.0759 


32 


b-W 


1750 


41.714 (0.132) 


0.290301 (21) 


0.4 


0.0759 


44 


W 


60 


41.948 (0.927) 


0.290283 (52) 


0.35 


0.075628 


48 


w 


130 


58.699 (0.420) 


0.255800 (18) 


0.3 


0.0754 


32* 


S-W 


345 


87.449 (0.758) 


0.220540 (75) 


0.3 


0.0754 


48 


w 


406 


87.821 (0.555) 


0.220482 (19) 


0.275 


0.075313 


48 


W 


53 


104.156 (1.305) 


0.204771 (40) 


0.25 


0.075231 


60* 


w 


42 


130.798 (1.369) 


0.188119 (31) 


0.2 


0.0751 


48* 


w 


27 


203.828 (3.058) 


0.156649 (103) 


0.2 


0.0751 


52^ 


w 


48 


201.191 (6.140) 


0.156535 (65) 


0.2 


0.0751 


60^ 


w 


7 


202.398 (8.614) 


0.156476 (148) 


0.15 


0.074968 


68^ 


w 


25 


460.199 (4.884) 


0.112611 (51) 


0.1 


0.0749 


68^ 


w 


24 


1125.444 (36.365) 


0.077358 (123) 


0.1 


0.0749 


72^ 


w 


8 


1140.880 (39.025) 


0.077515 (210) 



Table 1: The details of the lattice simulations for each k corresponding to minput- In the algo- 
rithm column, 'S-W stands for the Swendsen-Wang algorithm [15], while 'W stands for the Wolff 
algorithm [16]. 'Ksweeps' stands for sweeps multiplied by 10'^. 



shall illustrate in the following, their conclusion does not apply: to check the true behaviour 
of = 2Km'jjXia.tt one should first identify correctly the operative definition of mn on the 
lattice. The mass values reported by Balog et al. are not reliable determinations of the 
physical Higgs mass if one requires the theoretical consistency of the adopted definition of 
'mass'. To fully appreciate what is going on, we need to go back to Refs. [5-7] (otherwise 
ignored by the authors of Ref . [8] ) . 

In those calculations, one was fitting the lattice data for the connected propagator to 
the (lattice version of the) two-parameter form 

GM = ^^P^- (4.3) 

This is a clean and simple strategy: there is no reason to restrict the analysis of the 
propagator to the limit p — > 0. In fact, in a free-field theory (the lattice version of) 
Eq.([4.3|) is valid in the full range < < (with A vr/a). In a 'trivial' theory 
one expects a two-parameter fit to the propagator data to have small residual corrections. 
These, however, should become smaller and smaller by approaching the continuum limit. 

In this way, after computing the lattice zero- momentum susceptibility xiattj it becomes 
possible to compare the measured value of = 2Kmj^3^^^xiatt with the fitted Zp^op, both in 
the symmetric and broken phases. While no difference was found in the symmetric phase 
(see Fig. 3), Z^ and Zprop were found to be sizeably different in the broken phase. The 
discrepancy was found to become larger and larger by approaching the continuum limit 
and thus cannot be explained in terms of residual perturbative corrections to Eq.( [4.3D . In 
particular, Zprop was very slowly varying and steadily approaching unity from below in the 
continuum limit consistently with Kallen-Lehmann representation and 'triviality'. Z^p, on 
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the other hand, was found to rapidly increase above unity in the same hmit. The observed 
trend was consistent with the logarithmically increasing trend predicted in Eq.(2.9). 

This conclusion was based on the values of miatt extracted by skipping the lowest 
3-4 p^— values in the fit to the propagator data. In fact, differently from the symmetric- 
phase simulations, where Eq.([4.3|) reproduces the data in the full range 
in the broken-symmetry phase the two-parameter form Eq. ( [4. 3D does not reproduce the 
propagator data down to p = 0, see Fig. 4 (and Figs. 3, 5, 6 of Ref. [6]). Thus one has to 
choose. Either a) to restrict to the lowest 3-4 p"^— data, and obtain from the fit a pair of 
values {Ziov,, miow) or b) skip these first few data and fit the much larger sample of higher 
momentum data thus obtaining the other pair (.^highi "^high)- 

The two sets of mass values differ sizeably. For instance for k = 0.076, 0.07512, 0.07504 
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Figure 3: The lattice data for the re-scaled propagator in the symmetric phase at k = 0.074 as 
a fmiction of the lattice momentum = 2sinp^/2. The zero-momentum full point is defined as 
— {2 K)mf ^^^Xintt- The dashed line indicates the value of .^prop- Data are taken from Ref. [6]. 



the results of Refs. [6] were respectively mtigh = 0.4286(46), 0.2062(41), 0.1723(34). On 
the other hand, the alternative values from the lowest momentum data are respectively 
miow = 0.392(4), 0.1737(24), 0.1419(17) for the same k's. 

Numerically, the typical miow's are extremely close to the other mass definition mTs(k = 0), 
the mass extracted at zero 3-momentum from the exponential decay (TS='Time Slice') of 
the connected two-point correlator 

Ci(t,0;k) = (5,(t;k)5e(0;k) + 5,(t;k)5,(0;k))eonn, (4.4) 

where 

5e(t;k) = J];0(x,t)cos(k-x), (4.5) 

X 

5,(t;k) = J^0(x,t)sin(k-x). (4.6) 

X 

Here, t is the Euclidean time; x is the spatial part of the site 4- vector x^; k is the lattice 
momentum k = (27r/L)(na;, n^), with {nx,ny,nz) non-negative integers; and (...)conn 
denotes the connected expectation value with respect to the lattice action, Eq. (|3.1| ). In 
this way, parameterizing the correlator Ci in terms of the energy as {Lt being the lattice 
size in time direction) 



Ci(t,0;k) = A[eM-Ekt) + exp{-Ek{Lt -t))], 



(4.7) 
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K=0.07512 



_6 9. — ^ a s — s ^ 



"■'0 1 2 3 4 5 6 7 8 

A2 
P 

Figure 4: The lattice data for tiie re-scaled propagator in the broken phase at k = 0.07512 as 
a function of the lattice momentum = 2sinp^/2. The zero- momentum full point is defined as 
Zip = (2K)rn^j^j^xiatt- The dashed line indicates the value of Zprop- Data are taken from Ref. [6]. 

the mass can be determined through the lattice dispersion relation 

3 

m^s(k) = 2(cosh£;fc - 1) - 2 ^ (1 - cos A;^) . (4.8) 

In a free-field theory mxs is independent of k and coincides with m-iatt from Eq. ( ^.31 ) . 

Now, using the corresponding susceptibility values xiatt =37.85(6), 193.1(1.7), 293.4(2.9) 
reported in Ref. [6] and the above values of niH = m^igh one obtains a set of logarithmically 
increasing values (see Fig. 5): 2,^ = 

(2K)m^Xiatt = 1.05(2), 1.23(5), 1.31(5) as predicted in 

Refs. [2-4]. 

This is in contrast with the values of the other quantity Ziow ~ = {2K)m'^^^xia.tt = 
0.884(18), 0.875(24), 0.887(21) which remain remarkably stable (see the corresponding en- 
tries for Zn shown in Table 4 of Ref. [8] ) . 

Therefore, the crucial question raised by the propagator data is the following. Which 
is the 'true' lattice definition of niH ? Is the p ^ choice (adopted by Balog et al.) and 
represented by miow ~ w.Ts(k = 0), or that obtained from mhigh, as proposed in Refs. [6]? 
As discussed in Ref. [6], there are several arguments that suggest the correct definition of 
rriH to be obtained from mhigh; thus regarding the other values extracted from the very 
low-momentum region as a symptom of the distinct dynamics of the scalar condensate. 
Here we list a few: 

i) in the continuum theory, the shifted fluctuation field is defined as the 
7^ projection of the full quantum fleld. However in a lattice simulation with periodic 
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Figure 5: The two rescaling factors Zp^op and obtained in Ref. [6] as a function of miatt- 



boundary conditions, where the momenta pi are proportional to an integer number rij 
times the inverse lattice size 1/L, the notion 7^ is ambiguous. In fact, any finite 
set of integers will evolve onto the = state in the limit L — > 00. This means that, 
for a given lattice mass, by increasing the lattice size, to separate unambiguously the 
genuine finite-momentum fluctuation field from the 'condensate' itself, one should increase 
correspondingly the set of integers nj. This gives a clean physical meaning to the fit 
obtained by skipping the lowest momentum propagator data and to the pair of parameters 

ii) the values obtained in Ref. [6] (respectively for k =0.076, 0.07512,0.07504) 
"ihigh =0.4286(46), 0.2062(41), 0.1723(34) give a physical mass mn that scales as expected. 
This can easily be checked from the remarkable agreement between these values and the 
predicted trend Eq. (|4.1| ) which is valid at the leading- log level both in perturbation theory 
and in the alternative picture of Refs. [2-4], see Eqs.( |2.1^ ) and ( p. 15 ). In this way one 
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Figure 6: The time-slice mass mxs (k) at different values of the three- momentum in the symmetric 
phase. The shaded area represents the value mutt = 0.2141(28) obtained from the fit to the 
propagator data in Fig. 3. The black square is the result of Ref. [19]. 



finds A = 17.282(364), Kc = 0.074836(15) and predicts nin ~0.499, 0.408, 0.293, 0.199 
for K =0.0764, 0.0759, 0.0754, 0.0751 in good agreement with the values 0.5, 0.4, 0.3, 0.2 
reported in the LW Table. 

iii) the identification of miow ~ rriTsi^ = 0) as the physical ruH contradicts the 
observed dependence of mTs(k) on |k| in the limit k ^ 0. In fact, in the symmetric phase 
(see Fig. 6) the energy spectrum has the form Vk^ + m? up to momenta k^ ~ 25m^(!). 
However in the broken phase the energy is not reproduced by the (lattice version of the) 
form \/k^ + const (see Fig. 7) and thus the very notion of 'mass' becomes problematic. At 
larger k, where mTs(k) becomes insensitive to k, it agrees well with mhigh- The very 
different behaviour between symmetric and broken phase shown in Figs. 6 and 7 has no 
counterpart in the conventional picture. 

iv) the values of Zprop obtained from the fit to the higher-momentum data in 
Ref. [6] for k = 0.076,0.07512,0.07504, Zhigh =0.9321(44), 0.9551(21), 0.9566(13), exhibit 
a monotonical increase toward unity (from below) as expected on the base of the Kallen- 
Lehmann representation in a 'trivial' theory. This confirms that fitting the propagator 
skipping the lowest momentum data gives consistent results. On the contrary, the values 
for Zio„ ~ Zji obtained from the fit to the lowest momentum points alone, remain constant 
to ~ 0.88 in the limit k ^ Kc (see also the values in Table 4 of Ref. [8]). Assuming this 
definition as the correct one to be used in the Kallen-Lehmann representation, one would 
find an inner contradiction with the non-interacting nature of the shifted fluctuation field 
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Figure 7: The time-slice mass mxs (k) for several values of the spatial momentum in the broken 
phase. The zero- momentum value is our result 0.392(4) obtained from a simulation of 6 Msweeps 
on a 20^ lattice (in perfect agreement with the value 0.392(1) reported in Table 4 of Ref. [8]). 
The shaded area represents the value miatt — 0.42865(456) obtained from the fit to the higher 
momentum propagator data. 



in the continuum limit. 

No trace of this discussion is found in Ref. [8]. These authors, ignoring the evident 
difference between Fig. 3 and Fig. 4, as well as between Fig. 6 and Fig. 7, do not address 
the physical consistency of the mass definition extracted from the very low-momentum 
region of the lattice data. They just limit themselves to the remark that, for the A<I>^ case. 



the lattice propagator is consistently reproduced by Eq. (|4.3D for p — > 0. However, this 
agreement is not relevant here since the Ising limit is known to anticipate much better the 
true aspects of the theory in finite lattices. In the Ising case, in fact, the fit is not good 
and Balog et al. are forced to extract the lattice mass by fitting just to the three lowest 
momentum points, precisely as with the parameter miow of Ref. [6]. 

It is also surprising their claim (see the Abstract of Ref. [8]) that the lattice data are 
consistent with the perturbative predictions. We have seen in the previous section that, if 
one requires consistency between 2-loop predictions and the data for t^^X2(0), one has to 
replace the integration constant C'^ given in the LW's Table 1 with the value = 6.67(4) 
given m Eq.(p^). Therefore, using the perturbative relation reported in Eq.(2.14) of 
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Ref. [8], namely {or = t|^) 



ZR = i2K)C',{l--aR + 0{a%)) (4.9) 

one can compare this Zr with the other estimate extracted from the product m|jX2(0) and 
reported in Table 4 of Ref. [8]. 

In this case, for k = 0.076 (where the relevant value is qr = 30.37(28)) Eq.( |4.9| ) predicts 
Zr = 0.976(6) while the value reported in Table 4 by Balog et al. is Zr = 0.896(4). 
Analogously, for k = 0.0751 (where the relevant values is gR = 20.51(74)) Eq. (|4.9D predicts 
again Zr = 0.976(6) while the value reported in Table 4 by Balog et al. is Zr = 0.883(17). 
These discrepancies, that are at the level of ~ lOcr and ~ 5a, can hardly be considered 
indicative of theoretical consistency. 

Quite independently, the authors of Ref. [8] have not shown that the tur^s reported 
in their Table 4 lie on well defined RG trajectories niR = mR{K) using the same Kc = 
0.074833(17) obtained from the 2-loop fit to the data for f^X2(0). In particular, this 
means that their value niR = 0.395(1) for k = 0.076 has to come out consistent with the 
other value tur = 0.1688(15) for k = 0.0751. To clarify this issue, they should quote the 
new values of the mass that replace the LW entries for all values of k. 



5. Summary and conclusions 

Let us now try to summarize the various theoretical and numerical points addressed in 
this paper. We shall start by observing that the field rescaling is usually viewed as an 
'operatorial statement' between bare and renormalized fields operators of the type 

" $b(x) = ^$ij,(x) " (5.1) 

As pointed out in Ref. [12], this relation is a consistent short-hand notation in a theory 
where the field operator admits an asymptotic Fock representation, as in QED. In the 
presence of spontaneous symmetry breaking it has no rigorous basis since the Fock repre- 
sentation exists only for the shifted fluctuating field /i_b(x) = ^six) — {^b), the one with 
a vanishing expectation value. 

For this reason, following Refs. [2-4] (see the discussion given in Section 2), one can 
consider theoretical frameworks that are fully consistent with 'triviality' but where such 
an operatorial relation is not valid. In fact, the physical conditions used to determine the 



hs — hR relation, through du Z — Zpj-Qp ~ 1 in Eg . ( |274|) , can be basically different from those 
used in the vb — vr case, through a Z = ~ InA in Eq. (|2.6| ). In this case, if one wants 
to introduce a renormalized field operator ^r{x) (whose vacuum expectation value is what 
one defines by vr and whose fluctuating part is what one defines by Hr), this cannot be 
related to $s(x) by means of Eq. (|5.lD . 

As reviewed at the end of Section 2, this 'subtlety' has a substantial phenomenological 
implication: 'triviality', by itself, cannot be used to place upper bounds on the Higgs boson 
mass. Thus, the importance of the issue requires dedicated lattice simulations to check the 
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validity of the logarithmic trend predicted in Eq.(2^) by measuring the relative scaling of 
the physical Higgs mass squared m'jj vs. the zero-momentum susceptibility. 

Clearly, the answer to this question depends on the given procedure adopted on the 
lattice to extract mn = ?7iiatt- Our point is that a correct choice requires to fulfill several 
consistency tests, taking into account the free- field nature of the fluctuation field in the 
continuum limit. Thus one should first find a 4- momentum region where the connected 
propagator is described by the (lattice version of the) two-parameter form 



One should also check that the fitted value mn = ?Tiiatt controls the energy eigenvalues 
E = E{k) governing the exponential time decay of the connected correlator in a given 
region of 3-momentum k. In fact, where the energy is not reproduced by the (lattice 
version of the) form Vk^ + const., the very notion of mass becomes problematic. Finally, 
one should also check that the fitted values of Z^^-op approach unity from below in the 
continuum limit as required by the Kallen-Lehmann representation in a 'trivial' theory. 

Now, the usual assumption is to extract mu = fn\s,tt from the p ^ Q and/or k ^ 
limits. This is certainly valid in a simulation performed in the symmetric phase (see Figs. 3 
and 6) where the whole 3- and 4- momentum regions give the same indications. However, 
in the broken-symmetry phase, where the vacuum is some sort of 'condensate', there might 
be reasons to avoid the strict zero-momentum limit to extract the physical particle mass. 

For instance, as mentioned in Sect. 4, in the continuum theory, the shifted fluctuation 
field is defined as the Pp, ^ projection of the full quantum field. However in a lattice 
simulation with periodic boundary conditions, where the momenta are proportional to an 
integer number rii times the inverse lattice size 1/L, the notion 7^ is ambiguous. In fact, 
any finite set of integers will evolve onto the p^ = state in the limit L — > 00. Therefore, 
for a given lattice mass, by increasing the lattice size, to separate unambiguously the 
genuine finite-momentum fluctuation field from the 'condensate' itself, one should increase 
correspondingly the set of integers rij. 

In addition, there are precise physical motivations suggested by the non-relativistic 
limit of a broken-symmetry theory: the low-temperature phase of a hard-sphere Bose 
gas [20]. In this case, the low-lying excitations for k — > are phonons, i.e. collective 
oscillations of the hard-sphere system whose energy grows linearly £'ph(k) ~ Cs|k|, 
being the speed of sound. Only at larger |k| does the energy spectrum grow quadratically. 
Therefore, a determination of the effective hard-sphere mass through the non-relativistic 
1-particle relation £'i_part(k) ~ 2m g cannot be obtained from the k ^ limit of the energy 
spectrum which is dominated by the phonon branch. 

This type of problems were preliminarily considered in Ref. [6]. The result of that 
investigation was that the physical mass mn = miatt is obtained from the propagator data 
after skipping the lowest 3-4 momentum points. The mass value m-iow) obtained from the 
lowest momentum points, which is numerically close to the other mass extracted from the 
exponential decay of the connected correlator at zero 3-momentum, does not fulfill the 
same consistency checks. 
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Now, for K = 0.076, 0.07512, 0.07504 the results of Ref. [6] were respectively mn = 
miatt = 0.4286(46), 0.2062(41), 0.1723(34). In this way, using the susceptibility values 
Xi^tt =37.85(6), 193.1(1.7), 293.4(2.9) reported in Ref. [6] one obtains a set of logarith- 
mically increasing values: = (2/t)m|^xiatt = 1.05(2), 1.23(5), 1.31(5) as predicted in 
Refs. [2-4] (see Fig. 5). This is in contrast with the values obtained using m-iow for which 
the quantity Zr = {2 K)mf ^^xia.tt ~ 0.88 remains remarkably stable. 

To provide further evidence, we replaced in Ref. [1] the direct evaluation of niH = ?7iiatt 
with the theoretical input values predicted in the LW Tables. This is not in contradic- 
tion with the previous strategy, in fact the values obtained in Ref. [6] miatt =0.4286(46), 
0.2062(41), 0.1723(34) give a physical mass rriH that scales as in Eq.(^T]) and that is in 
good agreement with the values tuh = ?7iinput(^) reported in the LW Table. Using the 
values of xiatt reported in our Table 1 and the LW entries for the mass, this additional test 
confirms that the quantity Z^p = 2Kmf^p^^xia.tt increases logarithmically when approaching 
the continuum limit (see Fig. 2). 

Now, Balog et al., being aware that the above numerical evidences have "..serious non 
standard implications for the Higgs sector of the Standard Model" (see the Conclusions of 
Ref. [8]), have performed a new analysis. They claim that the LW entries for m/j, should 
be replaced by new values (whose consistency with RG trajectories tur = mji{K), however, 
has not been shown). As far as we can see, they have essentially re-discovered the result 
of Ref. [6] that the quantity Zr = {2K)mf^^xia.tt is a constant ~ 0.88. However, ignoring 
Ref. [6], they fail to appreciate why their values tur = miow do not represent a consistent 
lattice definition of mn- 

A point where we accept their criticism concerns the lattice data for f^X2(0). Restrict- 
ing to this observable, a definitive test of the leading-logarithmic trend has to be postponed 
to data taken closer to the critical point where the non-leading terms associated with the 
scale of the logs should become unessential. 

However, as discussed at the end of Sect. 3, within the perturbative framework, the 
data for t'^X2(0) give precise information on the integration constants C- (i=l,2,3) that 
should be used for the matching with the symmetric phase. Using the precise outcome of 



the fit C2 = 6.67(4) Eq. (O) predicts Zr ~ 0.976(6) with 5-lOa discrepancies with respect 



to the values reported by Balog et al in their Table 4. 

This is precisely the same discrepancy pointed out by Jansen et al. Ref. [21] whose 
origin cannot be understood ignoring the main point of Ref. [6] and of this paper. In the 
broken phase, a naive zero 4- momentum limit of the connected propagator and/or a naive 
zero 3-momentum limit of the energy eigenvalue controlling the exponential decay of the 
connected correlator, do not provide consistent estimates of the physical mass parameter 
rriH- The phenomenological implications are substantial. Once the quadratic shape of the 
effective potential (the inverse susceptibility) , that is a pure zero- momentum quantity, does 
not scale uniformly with rnj^, although the finite- momentum fiuctuation field becomes free 
in the continuum limit, one cannot use 'triviality' to place upper bounds on mn- 
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